Two-orbital Schwinger Boson Representation of Spin-One: Application to a 
Non-abelian Spin Liquid with Quaternion Gauge Field 
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A non-abelian spin liquid in triangular lattice spin-1 systems was recently formulated in the form 
of continuum field th eory [T. Grover, a nd T. Senthil, Phys. Rev. Lett. 107, 077203 (2011); Cenke 
Xu, A.W.W. Ludwig, arXiv: 1012.5671 . It has spin-1/2 bosonic spinons coupled to emergent quater- 
nion gauge fields, and can be obtained by quantum disordering a non-collinear spin nematic order 
hypothesized to describe NiGa2S4 [H. Tsunetsugu, and M. Arikawa, J. Phys. Soc. Jpn. 75, 083701 
(2006)], However a microscopic lattice description, e.g. the lattice spinon (mean-field) Hamiltonian 
and the spin wavefunction, has been missing, and it has been noted that the standard Schwinger 
boson or bosonic triplon representations of spin-1 cannot describe this spin liquid. In this paper a 
two-orbital Schwinger boson representation for spin-1 systems is developed and used to construct 
a mean-field description of this quaternion spin liquid. Projecting the mean-field state produces a 
prototype wavefunction, which is a superposition of close-packed AKLT loop configurations with 
nontrivial amplitudes. This new formalism and related wavefunctions may be generalized to higher 
spin systems and can possibly produce spin liquid states with even richer emergent gauge structures. 
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References llO 

Spin liquid states in more than one spatial dimensions 
were proposed more than three decades agoi. They are 
ground states of Mott insulators with no spontaneous 
symmetry breaking, thus beyond the symmetry break- 
ing paradigm of phases. Many parent Hamiltonians of 
spin liquids wither— and without 7 ^ spin SU (2) symmetry 
have been constructed. Extensive numerical studies on 
semi-realistic models have shown evidences of spin liquid 
ground states in quantum spin models on triangular— ~— 
and kagome lattices 1 ^ - — and also in electronic Hubbard 
model s 17 ' 18 . In the last decade several promising candi- 
date materials have also emerged, a review of which is 
given in Ref. [l9|. 

One way to understand some of the spin liquid states 
is by disordering a spin SU (2) symmetry breaking order 
without proliferating topological defects 2 ^— . Low en- 
ergy theory of such description usually contains gapped 
bosonic spinon and emergent gauge field, and the phase 
transition from quantum disordered (spin liquid) state 



to ordered state is the condensation of the bosonic 
spinon22r— . This approach is believed to work better in 
the deep Mott insulating limit, where a short-range quan- 
tum spin model is appropriate. To get quantum spin liq- 
uid ground state, the conventional wisdom suggests that 
low spin value is important, and spin-1/2 is the best. 
For spin-1/2 systems the only single-site 5'C/(2)-breaking 
order parameter is the local magnetic dipole vector (S) 
(which also breaks time- reversal symmetry). Long-range 
order of (S) is the most commonly used starting point of 
this bosonic spin liquid approach 2 ^—. 

However spin-1 systems may have magnetic 
quadrupole order (spin nematic order hereafter) 
that breaks spin SU (2) symmetry with zero local dipole 
moments and no time-reserval symmetry breaking. The 
spin nematic order parameter is the real symmetric 
traceless matrix Q ab = {{S a S b + S b S a )/2) - (2/3)5 ab 
(a, 6 = x,y,z). In this paper only the uniaxial spin 
nematic order will be considered, which can be described 



by the "director" n as Q 



ab 



(1/3)6" 



Spin nematic orders have been propose d 24 ' 25 for the 
spin-1 triangular lattice material NiGa2S4 which had 
some experimental evidence of a ground state without 
magnetic dipole order—. In particular Tsunetsugu and 
Arikawa24 proposed an interesting three-sublattice spin 
nematic order, with the directors on the three sublat- 
tices perpendicular to each other (see Fig. [T|) . This state 
was also found in a numerical study of spin-1 nearest- 
neighbor bilinear-biquadratic Hciscnberg model on tri- 



angular lattici 



.27 



Very recently two group 



,28.2!) 



sidered possible spin liquid states by disordering this 
"Tsunetsugu- Arikawa state" ( "antiferroquadrupolar or- 
der" in Ref. [27l "tetrad order" in Ref. |29[ ) , and found an 
interesting non-abelian spin liquid with spin-1/2 bosonic 
spinons coupled to emergent gauge fields in the quater- 
nion group Qg, a discrete non-abelian group with eight 
elements defined as Qs = {±1, ±i, ±j, ±k} with the mul- 
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FIG. 1: (Color online) Triangular lattice with the three sub- 
lattices labelled as A, B, C respectively. 61,2,3 are three lattice 
translation vectors. <j s and ad illustrate the two reflection gen- 
erators of point group Ce v ■ An example of the three-sublattice 
spin nematic ordep^ is illustrated, with the directors x, y, z 
drawn as red horizontal bars, green vertical bars, and blue 
dots, respectively. The quaternion IGG elements on A, B, C 
sites are the A, B, C representations of Qg in Eqs. ()25aH25c]) 
respectively. Arrow from site r to r' means the mean-field 
ansatz A rr > = -Aj,, r = U — (1 - ir x - ir v + ir z )/2. 

tiplication rule i 2 = j 2 = k 2 = ijk = — 1. The continuum 
field theory for this spin liquid and for the ordering tran- 
sition, and many interesting properties, e.g. topological 
degeneracy, have been worked out in Refs. I28ll29l How- 
ever no microscopic lattice spinon (mean-field) Hamilto- 
nian was given and it was not clear how to construct a 
spin wavefunction for this "quaternion spin liquid" . The 
authors of Ref. [H realized that it is impossible to de- 
scribe this spin liquid by the Schwinger boson or "bosonic 
triplon" construction, and suggested that "such a spin 
liquid cannot be obtained by the standard projective con- 
struction for spin liquids"^,. In this paper a new and non- 
standard projective construction in terms of a two-orbital 
Schwinger boson representation of spin-1 will be formu- 
lated, and used to construct a mean-field description of 
the quaternion spin liquid and its ordering transtion to 
the Tsunctsugu-Arikawa state, and produce a prototype 
spin wavefunction by projecting the mean-field state. 

Before proceeding to the new formalism it is worth re- 
viewing how the old Schwinger boson and bosonic triplon 
constructions fail. The Schwinger boson representation 
for a spin-S 1 system 3 ^ is to write the spin operators as (er 
are the spin Pauli matrices) 

s,s'=t4 

with the constraint ^2, s b\b s ~ 2S. This representa- 
tion has a U(l) redundancy, b —> e l6 b (i = 9 is 
real). For the gapped spin liquid phase to be stable in 
two-spatial-dimcnsion(2D) it is necessary to "Higgs" this 
U(l) gauge field to a discrete subgroup, usually Z2 by bo- 
son pairin g 20 ' 21 ' 31 . However the non-abelian quaternion 
gauge group is not a subgroup of this abclian U(l) gauge 
structure. Moreover it is impossible to describe the spin 
nematic phase in this formulation as a boson condensate, 
because a nonzero condensate z = (b) necessarily pro- 
duces a nonzero magnetic dipole moments M = z' ■ a ■ z jl 
with size \M\ = \z\ 2 /2. 



The bosonic triplon representatio n 32 ' 33 uses a three- 
component boson b and writes the spin operator as 

S = -itf x b, (2) 

with the constraint b' ■ b = 1. Both magnetic order and 
spin nematic order can be described by a boson conden- 
sate in this construction. However the gauge freedom 
is again U(l), b — > e b, with no quaternion subgroup. 
Moreover b is not a spin- 1/2 representation of spin SU (2) 
symmetry. What is needed for the quaternion spin liquid 
is a new representation of spin-1 by spin-1/2 bosons(like 
the Schwinger boson) with large enough gauge freedom, 
which can describe the spin nematic order semiclassically 
by boson condcnsation(likc the bosonic triplon). 

The outline of this paper is as follows. In section Q] 
the two-orbital Schwinger boson representation for spin- 
1 will be introduced. Generic mean-field theories of spin 
liquids in this representation and related gauge structure 
and generic (projected) mean-field wavefunctions will be 
presented. In section [TT1 a mean-field Hamiltonian for the 
quaternion spin liquid on triangular lattice will be con- 
structed and analyzed. The Tsunctsugu-Arikawa state 
will be obtained by boson condensation. A prototype 
spin wavefunction for the spin liquid will be presented. 
In section IIIII remaining issues and possible extensions 
will be discussed. 



I. TWO-ORBITAL SCHWINGER BOSON 
REPRESENTATION OF SPIN-ONE AND 
MEAN-FIELD THEORIES FOR SPIN LIQUIDS 

Spin-1 can be viewed as a symmetric combination of 
two spin-1/2. The Affleck-Kennedy-Licb-Tasaki(AKLT) 
state was originally defined in this way 3 -^ Use this 
"two-orbital" picture and introduce Schwinger bosons for 
each spin-1/2, the spin-1 operators in terms of the four- 
component bosons are 

2 2 

S = Y. S « = Y.2 £ ^ss'Ksn (3) 

a=l a=l s,s'="\-,l 

where a = 1,2 labels orbital, s,s' =t,| label spin. This 
bosonic representation for SU(N) spins were briefly men- 
tioned in Ref. From this alone the gauge freedom 
seems to be orbital U(2), namely b as — > J2p u apbp s with 

^2p w q/3 u q '^ = $aa'- However the gauge transformations 
should also leave the constraints invariant. 

The first constraint requires that the total number of 
bosons is two, 

n = J2n a =J2 b LKs= 2 - ( 4 ) 

a as 

Define the orbital pseudo-spins T as 

T = E \ b Lr ap b psl (5) 

8,<X,/3 
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where r arc orbital Pauli matrices. There are ten states 
for two bosons, one S = 0, T = spin-orbital singlet 
j2-3/2^f , T v a y . 6+) |0) where |0) is the boson vacuum] 
and nine S = 1,T = 1 states. The second constraint is 
to project out the singlet state, 



(6 f • r y a y ■ tf){b-T y a v -b) = 0. 



(6) 



The S = 1, T = 1 states can be arranged into a matrix, 
with row and column indices S z ,T Z = 1, 0, —1, 



V2 
\ V2 



\ 



(4t) : 

V2 I 



|0>. 



(7) 



The physical spin-1 states should be a linear combination 
of the columns in Eq. ([7]). This final constraint can be 
formally represented by 



Af-T = 0, 



(8) 



with complex unit vector Af (Af* ■ Af = 1). The 
chosen linear combination has coefficients — (Af x — 
iN y )/V2,N z , (Af x + iN v )/V2 respectively for the three 
columns in Eq. ([7]) . An equivalent form of Eq. ((SJ is 



E 

a,b—x,y,z 



T a ■ N a *N b ■ T b = 0. 



(9) 



The gauge freedom should leave this "vector Higgs con- 
densate" Af invariant up to a complex phase, or cquiva- 
lently leave the hermitian matrix Af a *Af b invariant. 

This parton construction unifies the conventional 
Schwinger boson Eq. (fTJ) and the two-orbital AKLT rep- 
resentation. For example, AT = (x + iy)/y/2 chooses the 
first column of Eq. (J7J) and is the old Schwinger boson 
representation. In contrast Af = z chooses the middle 
column of Eq. (J7J and is the AKLT representation. Three 
different cases of Af are discussed in the following. 

Case 1): Af* X Af = 0, then a complex phase can be 
chosen so that Af is a real vector. This represents an 
orbital nematic ordered state with the director Af, and 
is equivalent to the AKLT representation. The gauge 
freedom is U(l) x U(l) x Z 2 , generated by 



j»b, b 



3 #(a/"-t; 



b, b -> (Afx ■ t) • b, (10) 



where Afx is a real unit vector perpendicular to Af. Note 
that the last Z 2 does not commute with the second U(l) 
group, therefore scmidircct-product x is used. 

Case 2): Af* X Af ^ and Af ■ Af = 0, then the 
real and imaginary parts of Af are perpendicular to each 
other and of equal length y/l/2. This represents an or- 
bital dipole ordered state with the orbital moment along 
T = —iAf* x Af, and is equivalent to the single-orbital 
Schwinger boson representation Eq. (fTJ). The gauge free- 
dom is U(l), 



where 1 is the identity matrix. 

Case 3): Af* xAf^O and Af-Af^0, then a complex 
phase can be chosen so that the real and imaginary parts 
of Af are perpendicular to each other but of different 
length. The orbital dipole moment is nonzero along T = 
—iAf* x Af/\Af* x Af\. The gauge freedom is U(l) x Z 2 , 
generated by 



(T ■ r) ■ b. 



(12) 



Consider the Hciscnbcrg antifcrromagnetic interaction 
between two spin-1 at positions r and r' , namely S r ■ 
Sr' = E Q b Sra ■ S r 'fj- The right-hand-side can be 
Hubbard-Stratonovich decoupled in the same way as the 
Schwinger boson mean-field theory^. The mean-field 
Hamiltonian contains spin singlet boson pairing terms, 

Ar'.aM^A'M ~ ^A',^ + k - C -> OT 111 ShOTt foml 

b\ ■ A rr , (X) i(J v ■ b\, + h.c. where A rr ' is a generic 2x2 ma- 
trix in the orbital space. Note that A r , r = —Aj r , where 
superscript T stands for matrix transpose. For simplicity 
the boson hopping terms ^ s b\, as b r , a s are ignored. 

The constraints must be included in the mean-field 
theory by introducing onsite terms with Langrange 
multipliers^. The first constraint Eq. Q can be incorpo- 
rated by a chemical potential /i as fj,(n — 2). The second 
one Eq. © may be included as - A(6 f • T y a y ■ )(b ■ r y a y ■ b) 
with real A, however this is not quadratic in terms of 
bosons. A non-rigorous Hubbard-Stratonovich decou- 
pling can be performed to reduce this term to —r)(W ■ 
T y a v . &t) - v *(b- r y a y ■ b) + M 2 /A with a complex Lan- 
grange multiplier 77. This procedure may be made rigor- 
ous regardless of the sign of A by the tricks of Ref . [3?]. The 
final constraint Eq. © may be included as —^'Af-T+h.c. 
with a complex Lagrange multiplier p! . In summary the 
generic mean-field Hamiltonian for spin liquids is 

#mf = - ^2\bl ■ A rr , <X> ia y ■ b\, + h.c] + ^ fi r n r 

r,r' r 

-YJWM-T v (7 y -bl) + h. C ] 
r 

- y^[&^ • K(/x^AT r ) • r • b r ] + const., 

r 

(13) 

where 5ft means real part. The mean-field constraints are 
(n r ) MF =2, (&t . T V • &t) MF = 0, 



(bi-(Af r -T)-b r ) MF =0. 



(14) 



b _> e mi+*-r)h 



(11) 



The mean-field Hamiltonian Eq. (|T3|) is not gauge 
invariant. Under a site-dependent gauge transfor- 
mation b r — > G(r) ■ b r , the mean-field ansatz 
{p r , T) r , p' r , Af r , A rr '} should transform as 

Arr' G(r) ■ Arr' ■ G T (r'), r) r ->■ e 2W Vr , 

Hr + $t(fJ, r Afr) ■ r -> G(r) ■ [ Mr + %t{p r Afr) ■ t] • G-\r), 

(15) 
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where 9 defines the orbital-independent U(l) subgroup 
in Eqs. (|10II12[) . Gauge- invariant fluxes can be defined in 
analogy to the Schwinger boson or Sp{N) boson theory 3 ^. 
The loop expansion for the mean-field ground state en- 
ergy can also be performed, and a "flux expulsion" argu- 
ment for Hciscnberg models may be raised as well^. 

A realistic spin-1 Hamiltonian may contain the bi- 
quadratic interactions (Sy- S r >) 2 and multiple-spin inter- 
actions, which cannot be simply decoupled into quadratic 
terms of bosons. In this situation it is better to view 
the mean-field theory as a variational approach. The 
mean-field ground state after projection to physical spin- 
1 space can be used as a variational wavefunction. This 
viewpoint will be adopted throughout this paper, so no 
self-consistent equation of A will be solved, and the over- 
all scale of the ansatz does not matter. 

The mean-field ground state |MF) is generically 



|MF)=exp [\^l-fr 



T V -bl 



|0), (16) 



where f rr i are 2x2 matrices in the orbital space, and 
have the same symmetry and gauge transformation rule 
as the mean-field ansatz A rr ', e.g. f r , r = —f^ r ,- r = 



r' term with f rr = -f} 



r y is allowed but creates 



only onsite spin singlet and will be projected out. Each 
frr', a f}( b at b li _ b ai b lt^ tcrm creates a spin singlet from 
two spin-1/2 on bond rr' . The projection onto physical 
spin-1 space requires two bosons on every site and the 
onsite symmetrization of the two orbitals. The projected 
wave function V | MF) is therefore a superposition of close- 
packed (every site is covered once) loop configurations 
{£}, and on each loop £ the spin-1 form an AKLT state, 



V\MF) = Y[Wi\AKLT on I). 

{i} i 



(17) 



The "close-packed" loop configurations may involve 
bonds longer than nearest-neighbor. The amplitude fac- 
tor Wi for a length- L loop £ = (r\r 2 ■ ■ - r L ) is 



W t = (3/4) L / 2 • N L ■ Tr[r r J rir2 T r J r2r 



Jr 



(18) 



where Tr means matrix trace, Nl = yl + 3 • (—3) L , 
and T r = —ir v ■ (Af* ■ r) comes from the contraint Eq. (|8]) 
(r r = t x when Af r = z). The factor (3/4) L / 2 • N L is 
the overlap between one spin-1/2 dimer pattern and the 
AKLT stated. (3/4) i/2 produces an overall factor for 
the wavefunction and can be omitted. Nl ~ 1 when 
L is large. Wi is gauge invariant up to a global factor, 
due to the fact that G T (r) ■ (-ir v ) ■ (Af* ■ r) • G(r) = 
e 2W (-t T y) • G- X (r) • (Af* ■ t) ■ G(r) ~ (-»>») • (Af* r ■ r) 
up to a complex phase, for any G(r) in the gauge group 
Eqs. (|10H12|) . L can be 2 in which case the AKLT state 
is the "double-bond" spin singlet state formed by two 
spin-1. 



Wavef unctions for spin-1/2 spinon and gauge flux exci- 
tations can be constructed as well. The mean-field state 
with two spinons at r, r' is given by 

\B r , B r ,) M F = Tr[Bl ■ B]Tr[Bl ■ B r ,]\MF) (19) 
where B r y are 2x2 complex spinon state matrices, and 



= [ Jit J2 t 



(20) 



Projecting this state onto spin-1 space creates a super- 
position of configurations with one open AKLT chain 
l rr i from r to r' plus close-packed AKLT loops. The 
direction of the end-spins of open AKLT chain are given 
by Tt[B^<tB] (normalization requires |Tr[S^crS]| = 1). 
The amplitude for a length- (L + 2) open chain £ rr i = 
(rr\ . . . r^r') with the end-spins at r, r' in S z eigenstates 
s, s' =t, 4- is (L can be zero) 

W irr ,, ss > = (3/4) L / 2+1 • N L , SS , 

x (B r T r f rri T ri f rir2 . . . T rL f rLT ,T r iB r ,) ss ,, 

(21) 



witbM N L , a »> = V 1 + (-3)- L - 2 (a x - l) ss ,. Multiple 
(even) number of spinons can be constructed similarly. 
Gauge flux excitations and topological degeneracy of 
ground states will be demonstrated in the quantum limit 
of quaternion spin liquid in subsection III CI 

Spin liquids described in this way will have gapped 
spin-1/2 bosonic spinons. For them to be stable in 
2D it is necessary to "Higgs" the continuous com- 
pact gauge groups Eqs. (|10lll2p to a discrete subgroup. 
Many possibilities exist which can in principle be com- 
pletely classified by the projective symmetry group(PSG) 
metho d 36 ' 39 " — . This brute-force approach will not be 
attempted here, but the PSG language will be used to 
show that the mean-field theory indeed describes a spin 
liquid state with no symmetry breaking. This will be 
achieved by the explicit construction of the PSG el- 
ements, b r — > Gx(r) ■ b X { r ), for all generators X of 
the physical symmetry group (space group and time- 
reversal). The mean-field Hamiltonian shall be invariant 
under PSG actions. 

The case with uniform Af r = z will be considered here- 
after only, except subsection III Dl The physical spin-1 
states \S Z = +1, 0, —1) are 



\S Z = +1) = fclt^tio), 



(22a) 



1 

71 v 



0) = -^(6^6^ + 61^)10), (22b) 



IS 2 = -i) = 6iAi°>- ( 22c ) 

The constraints Eq. ((4]) and Eq. ([8]) can be rewritten as 
na = XXA s = 1 > a = 1,2. (23) 
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The U(l) x U(l) x Z 2 gauge group is, 

b -> e l9 (cos6 ■ 1 + sin (A • ir z ) • 6, or 

„ (24) 
6 e (cos</> • ir" 1 + sin0 ■ ir") • 6. 

For frustrated(e.g. triangular) lattices the orbital- 
independent U(l) freedom e will be removed by boson 
pairings in mean-field theory. The remaining U(l) xi Z 2 
[a subgroup of SU(2) by setting 9 = in Eq. (pi)) ] will 
be the starting point of discussions hereafter. It is non- 
abelian and contains the quaternion group. 

A semiclassical picture of the uniaxial spin nematic 
order is that the two spin- 1/2 have dipole moments 
antiparallcl to each other and along the direction of 
the director, (Si) = —(S2) oc n. This can be 
achieved by a single spin-orbital-cntangled condensate, 
e.g. ((&i t , 614., & 2 f, &24.)) oc (1,0,0,1) for n oc z. How- 
ever the quadrupole order parameter Q ab is naively zero 
because (S a ) is zero. This can be remedied by recogniz- 
ing that Q ab = (£a,*0SS^ + S b a S%))/2 - (2/3)S ab = 
(S^S^ + S b S%) - (l/6)S ab , and the last expression is 
nontrivial in this condensate state [although not in the 
form of n a n b — (l/3)S ab n 2 ]. This can be further justi- 
fied by projecting the coherent state from this conden- 
sate exp[w (b\^. + &2^)]|0) on to the physical spin-1 states 
Eqs. (|22aM22c[) . which gives the uniaxial spin nematic 
state (w 2 / \/2)\S z = 0) with director along z direction. 
In general the nematic director from a condensate (b) is 
given by (b)* ■ t z <t ■ (b). 

II. QUATERNION SPIN LIQUID ON 
TRIANGULAR LATTICE 

For a mean-field theory of quaternion spin liquid, the 
invariant gauge group(IGG) 36 i 39 must be a representa- 
tion of the quaternion group Qs = {±1, ±i, ±j, ±k}, 
with eight distinct IGG elements G q {r) £ U(l) xi Z 2 such 
that: 1) the ansatz {{i r , "Hr, ^, A/V = z,A rr '} are invari- 
ant under the actions Eq. (fT5j) of G q (r) for any q £ Qs, 
and no other element of Eq. (|2~i|) can do the same; 2) 
G q (r) is a representation of Qs for any site r. 

There are three distinct classes (labelled by A, B, C) 
of Qs representations on a single site. Representatives of 
each class are given below (G±i = ±1 for all classes), 

A: G ±i = Tir x , G±i = Tir y , G± k = T^ z ;(25a) 
B : G ±i = T*r\ G ±j = Tii~ z , G ±k = T^;(25b) 
C : G ±i = Tir z , G ±j = Tir x , G ±k = T »t*.(25c) 

Each class is generated by group conjugacy on its repre- 
sentative, G q -> G-Gq-G- 1 for G e t/(l)xZ 2 . Therefore 
by site(r)-dependent gauge transformations all G q (r) can 
be reduced to one of those in Eqs. (|25all25c|) . 

These Qs IGGs will constrain allowed ansatzs. For the 
onsite terms, the Qs IGGs demand fi' r = but put no 
constraint on r\ r . On the converse, however, fj,' r = and 
rj r do not reduce the U(l) x Z 2 freedom. 



Consider a bond rr' with nonzero A rr ' ■ The Qs IGGs 
demand A rr ' = G q (r) ■ A rr ' • G q (r') T for all q £ Qs- 
There are three possibilities for A rr ' depending on the 
representation choice combinations (rr'), 

( T y, (rr') = (AA),{BB),(CC), 
A^ oc I U, (rr') = (AB), (BC), (CA), (26) 
{ U T , (rr') = (AC),(BA),(CB), 

where the SU(2) matrix U = (1 - ir x - ir y + ir z )/2 
will appear frequently. Consider the converse problem, 
namely whether A rr ' can "Higgs" the gauge freedom 
to Qs- A rr ' oc t v will not do this job, because all 
G(r) = G(r') £ U(l) X Z 2 will keep A rr ' oc r* in- 
variant. The other two possibilities A rr > = U or U T 
will reduce the gauge freedom to Qs with the repre- 
sentation choices given above. For instance consider 
G(r') = cos0 • 1 + sincj) ■ ir z and A rr ' = U, the con- 
straint solves for G(r) = cos <fi ■ 1 + sin <j> ■ ir v which can 
be a member of the U(l) X Z 2 group only if is a in- 
tegral multiple of tt/2, restricting G(r) and G(r') to be 
members of Qs representations. 



A. Mean-field Theory on Triangular Lattice 

With the above general considerations a mean-field 
Hamiltonian of quaternion spin liquid can be constructed 
on the triangular lattice. Due to the three-sublattice 
structure it is natural to assign the three Qs represen- 
tations to the three corresponding sublattices. In this 
paper only the nearest-neighbor ansatz will be consid- 
ered, with A rr ' = U or —U T as shown in Fig. [1] Note 
that by the variational interpretation the overall scale 
of Arr' does not matter, and the overall complex phase 
of Arr' can be removed by a global orbital-independent 
U(l) phase rotation of bosons. Translation symmetry 
further requires uniform /i r = /! and r\ r = V- 

The mean- field Hamiltonian reads (up to a constant), 

Hmf = ^2 If* Hr ~ ' tV<jV ' b r + h.c.) 

r 

(27) 

- ^2(bt - U®i<7 v ■ bl +ea + h.c.) . 

d=l 

Physical symmetries are generated by two translations 
T\ y2 along ei t 2, two reflections a s and <Jd (see Fig. [TJ), 
and time-reversal T . T\ 2 and a s are trivial, ad reverses 
all bond orientations in Fig. [T] T changes the ansatz to 
their complex conjugate. Corresponding PSG elements 
are, 

Ti i2 : b r -> b r+ei 21 (28a) 
a s : b r ->• fe CTs ( r) , (28b) 

a d : b r ^y/T/2(ir x -iT y )-b ad(r) , (28c) 
T : b r ^ T y a y ■ b r . (28d) 
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The Q 8 IGG is defined in Eqs. (|25all25cj) for A, B, C sub- 
lattices respectively. Time-reversal symmetry restricts r\ 
to be real. With this PSG constructed the mean-field 
Hamiltonian describes a "symmetric spin liquid"— with 
no broken symmetry. 

The mean-field Hamiltonian can be solved in the 
same way as the Schwinger boson mean-field theories^!. 
Do the Fourier transform b k = N^J 2 J2 r e~ lkr b r 
where N sitc is the number of sites, and define Vffc = 

( fo fc,it' 6 fc,2T' 6 -fc.U' 6 -fc,2i) T ' Ec l- <E3> becomes 



Hmf = J2 H 



Hi -P k 

-pi mi 



(29) 



where P k = i(l - ir x + ir z ) 3£ fc - iT y ($t£ h +rj), 5? and 3 
are real and imaginary parts, = e lh ei + e lfc ' 62 
Do a singular value decomposition 



r,ike 3 



it 1 ;,)-'.' 



with U(2) matrices Uk,V k given by 
U k = 



(30) 



3+V3 



3-V^S 



c i7r/4 



3— v/3 



3— </3 



ITT I '4 



3+V3/ 



3+^ 



3+V3 



3-V3 



(31) 
(32) 



and real singular values pi, 2 (k) = V3^s£,k ± (5i£fc + V)- 
Define "Bogoliubov quasiparticlcs" 



( 7fe,it \ 
7*,2t 
7-fe,ii 



/Ci 5i \ 

c 2 o s 2 Uf/t 2x2 \ 

Si d (,02x2 F fe f J W< 

\o s 2 o cj 



(33) 

where Ci, 2 = V 1 + v/E h2 (k)/V2 and Si j2 = 
—pi,2(k)/Ei t 2 (fc)/2Ci,2, with the mean-field dispersions 



(34) 



Eq. ((29]) is diagonalized by this 51/(2,2) Bogoliubov 
transformation, 

#mf = [- B i( fe )(7l,it7fe,iT +7- fe ,ii7ife,u) 
fe 

+ ^2(fe)(7fe,2t7fe,2t + 7-fc,2;7lfc,2i) - 2 H • 

(35) 

The mean-field ground state energy per site is 

Emf = N~ it \ Y^[Ei(k)+E 2 (k)-2fi}. (36) 



The mean-field ground state is annihilated by all / 7fc jCKS , 
and is 



|MF) = exp 



[IE*-* 



y . ftt 



-fe 



|0) 



(37) 



where is the Fourier transform of f r 



Si/d 
S 2 /C 2 



h = -u k 



V, 



k ■ 



(38) 



and |0) is the vacuum of b bosons. 
The mean-field constraints arc 



WMF 



dE ] 



MP 



d 2 k[ 



2, (39) 



and 



0=(6 t -rV-i) t )HF = 



dE 



MF 



£l(fc) S 2 (fe) 

p 2 (fc) pi(k) 



dr) 



= / d 2 fc[^ 



-E 2 (k) E^k) 1 ' 

under the thermodynamic limit ^2 k — > J d 2 k 
I d 2 fc/(87r 2 /-\/3), where the integral is over the Brillouin 
zonc(BZ) with area 8n 2 /y/3. Note that the Qs IGG guar- 
antees that the boson numbers on the two orbitals are the 
same, so only the total density constraint is needed. 



B. Boson Condensation and Spin Nematic Order 

Similar to the standard Schwinger boson mean-field 
theorie a 20 ' 21 i 31 , the ordered state can be studied by re- 
laxing the total density constraint Eq. (|39[) . and driving 
the transition to ordered state by increasing boson den- 
sity. The minima of mean-field dispersions Eq. (|34[) arc 
always at BZ corners ±K, defined by K ■ e 2j 3 = —2n/3. 
So boson condensation will produce a three-sublattice or- 
der. 

Numerical solution of the mean-field critical point 
gives the critical ansatz parameters r\ c ps 0.5287 and 
fx c = 6 — rjc ~ 5.4713, the boson dispersions arc illus- 
trated in Fig. [21 The critical boson density n c sa 0.4492 
is however very low. Taken at face value this means the 
spin-1 system with n = 2 will be deep in the ordered 
phase. It is conceivable that the fluctuations ignored in 
the mean- field theory and/or farther neighbor couplings 
may stablize the spin liquid state. 

The boson condensate (VP) is a linear combination of 
the zero-energy eigenvectors of the mean-field Hamilto- 
nian Eq. (|29[) at the critical point. In real space it is 



1 



-3ir/4 p 7t/4. 



V4 C _1 



fcr/4 c 



e +iKr Wi 



w 2l 



(41) 
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/ \ M 

<^ r V - 



FIG. 2: Mean-field dispersions Eq. (|34[) along high symmetry 
lines at the critical point r\ c ~ 0.5287 and /i c = 6 — ry c ~ 
5.4713. Inset shows the Brillouin zone and high symmetry 
points and lines. Note that at K point there is a critical 
mode (Z in text) and a gapped mode (W in text). 



where w\ t 2 are complex coefficient, and constant c = 

y~2 + Note the eigenvectors at ±K (two rows in 
the first matrix) form a time-reversal pair. Define 



e 2iir/3 
e -vr/13 c 



eIT /12 c 
_ e -2i7r/3 



(42) 



and 5£/(2) rotor field Z 



'-2 



(43) 



The boson condensates (B r ) [B is defined in Eq. ([20]) ] on 
the three sublattices are 



B A ) = Z-(-ir v -Uf, (44a) 

B B ) = Z ■ (~ir y -U) T (-iT y -Uf, (44b) 
B c ) = Z. (44c) 



It is easy to see that the total dipole moment is 
(S) = (l/2)Tr[ZterZl] = 0. The nematic directors 
Ty[(B)^ct(B)t z ] on A, B 7 C sublattices are given respec- 
tively by, 



n A = -Tr[Z*o-ZT v ], n B = Tt[Z^(tZt x ], 
n c = Tr[Z^aZr z ], 



(45) 



which are perpendicular to each other. Therefore this is 
the Tsunetsugu-Arikawa spin nematic state. For exam- 
ple with Z = (1 — iT z )/y/2 the three directors are x, y, z 
on A, B, C sublattices respectively, which is the state de- 
picted in Fig. [I] The left SU(2) transformations of Z are 
spin rotations, the right SU(2) are sublattice-dependent 
orbital rotations [the gauge field is only f/(l)x[/(l)x Z 2 ]- 
The PSG transformation rules of Z can be derived from 



the PSG of lattice bosons Eqs. (|28aH28d|) . 



J 1.2 






Z 




z 


T 


z 


i 


z 


j 


z 


k 


z 



iT Z ) 



z\T 



1 



V2 

z\ 

Z-{-iT V f, 



(-IT V + IT Z ) 



z\T 



(46a) 

(46b) 

(46c) 

(46d) 
(46e) 
(46f) 
(46g) 



The symmetry allowed form of the low energy action 
would be (.D is covariant derivative), 



d 3 x(Tr[(D„Z^)(D v Z)] + m 2 Tr[Z f Z] 
+ uTr[Z t Z] 2 + ...}, 



(47) 



with 50(4) symmetr y 28 ' 29 . Other aspects of the field 
theory can be found in Refs. [28ll29l and will not bc re- 
peated here. 



C. Prototype Wavefunction in Quantum Limit 

The "quantum limit" is achieved by relaxing the total 
boson density constraint Eq. (j3"9")l and going to the low 
density limit 3 ^^ with (n) -C 1 and /j > 1. The mean- 
field constraint equations Eqs. (|39M0|) can be solved in 
power series of /a , 

1 35 

Mmf = 6m- 2 + — ^ 4 + 0(^ 6 ), 
r? = 6/j~ 2 + 81/i~ 4 + 0(fi- 6 ). 

By inverting the first equation every quantity can also be 
expressed in terms of (u)mf- 

The mean-field bond amplitudes f rr i in Eq. (|3"8"|) can 
also be expanded in power series of /i" 1 and will decay 
exponentially as /i~' r ~ r ' with respect to the distance |r— 
r'\. For example, bond amplitudes f rr i on the nearest- 
and second- and third-neighbor bonds are given by 



/r,r+ei 
fr ,r+2ei 



u 



rU 15/1 

r 3/i" 3 



-U 1 



4 

r 3/i~ 3 



16 



0(/x" 5 )], (49a) 



+ o( M - 5 )], 



(49b) 
(49c) 



and those related trivially by cyclic permutations of ei,2,3 
(three- fold rotations), and by f r , r = —,fj r ,- 

The wavefunction in the extreme quantum limit [i — > 
oo simplifies to the extreme "short-range resonating va- 
lence bonds" state, with nonzero amplitudes only on 
nearest-neighbor bonds. The overall factor of f rr > docs 



FIG. 3: Amplitudes [We in Eq. JTHJl without (3/4) i/2 ] of 
several short AKLT loops in the prototype wavefunction given 
by Eq. I|50p . The AKLT states reside counter-clockwise on the 
solid lines. Arrow from r to r' means the bond amplitude is 

frr' = ~fr'r = ~*W = ~i(t - ir* - ir y + i T Z )/2. 

not matter. The form of f rr i can be fixed by the Qg 
IGG without calculation and must be proportional to the 
mean- field ansatz Arr' — U or — U T . To simplify later 
discussions a factor —i is applied to f rr ', so 

{-ill, if r' = r + ei )2 ,3, 
iU T , if r' = r - ei j2 ,3, (50) 
0, otherwise, 

with U = (1 - it x - ir y + ir z )/2. This together with 
Eqs. (|17I18|) defines a prototype wavefunction which may 
describe the quaternion spin liquid (confinement is also 
possible). 

The loop weight We in Eq. (fP5|) becomes 

W e = N L • Tr[{-iT x )A ri r 2 . • . {-iT x )A rLri \. (51) 

with nearest- neighbor bonds < r\T2 >,...,< Tjj\ >, and 
(3/4) L / 2 factor ignored. The matrix product inside the 
trace symbol is a SU(2) matrix because every factor be- 
longs to SU(2), so the trace must be real. 

In fact the trace can only take three values, ±2 or 0. 
The proof is the following. Denote the number of bonds 
with orientation along the loop direction (Xnr i+ i = U) by 
JV+, and the number of those opposite to the loop direc- 
tion (Xnn + i = —U T ) by N-. Due to the three-sublatticc 
structure N + + 2N- = mod 3. This can be formally 
proved by assigning Z% numbers 0,1,2 to A, B,C sub- 
lattices respectively, and noting that travelling along(or 
against) a bond increase this number by unity(or two) 
modulo three. Use the identity —U T = {—it x )UU to re- 
place the N_ factors oiU T by 2iV_ factors of U, the ma- 
trix product becomes q\Uq<JJ. . . . qN + +2N-U, where the qs 
belong to the quaternion group {±1, ±ir x , ±ir v , ±ir z }. 
Use the commutation relations, U(±ir x ) = (^fiT z )U, 
U(±iT y ) = {±ir x )U, and U(±ir z ) = {Tir v )U, to move 
all the qi factors in front of all U factors, the matrix 
product becomes q ■ U N + +2N - where q is some quater- 
nion group element. Finally use U 3 = — 1 and the fact 
that iV + + 2iV_ is a multiple of three, the trace becomes 
Tr[±g] which can only be (if q is not ±1) or ±2. 

The AKLT normalization factor Nr, approaches unity 
in the long length L — > oo limit, so may be omitted with- 
out changing the long distance behavior especially the 



o 

FIG. 4: (Color online) 6x6 triangular lattice with periodic 
boundary condition. O indicates the origin site. The horizon- 
tal green (oblique blue) dash line, labelled by q~2 (<7i) and cut- 
ting through the top row (rightmost column) of green (blue) 
thick bonds, defines the creation operator of gauge flux in the 
non-contractible loop along e.2 (ei) direction. The red dash 
line cutting through several red thick bonds in the middle de- 
fines the creation operator of a pair of gauge fluxes at the two 
red dots. Gauge flux of class q £ Qs can be created by setting 
f rr i — > G q (r) ■ f rr ' for all bonds rr' cut by the defining dash 
line. The top-right cyan bond is shared by qi and 52- 



topological order. We for several short loops are pre- 
sented in Fig. [3J Interestingly the weight of "double- 
bond" singlet vanishes, so the wavefunction is made 
purely by AKLT loops of length L > 3. 

The prototype wavefunctions for gauge flux excitations 
and nontrivial topological sectors can be constructed in 
this quantum limit by the standard method (see e.g. 
Refs.|42j|43j). Gauge flux on a length-even loop (ri . . . T2l) 
can be denned as A rir2 {-A* r2r3 ) . . . A r2L _ ir2L {-A* r2Lri ). 
In the prototype wavefunction, the flux in each rhombus 
(unit cell) is 1. For a 6n x 6m lattice (m,n are integers) 
with periodic boundary condition the fluxes on the non- 
contractible (NC) loops are also 1. Creation operators of 
a pair of local fluxes are defined on the string on dual 
lattice connecting them. Creation operators of flux in 
a NC loop are denned on a NC loop of the dual lattice 
traversing it. Creation of gauge flux of class q e Qs 
amounts to A rr > — > G q (r) • A rr ' for all bonds rr' cut by 
the string or NC loop on dual lattice. Examples on a 6 x 6 
lattice are shown in Fig.fJ] Two fluxes on NC loops along 
(d = 1,2) direction can be explicitly defined as q& = 

•^0,5e d ( _ ^5e d) 4e d )-^4e d ,3e d ( _ -^3e £ i,2e £i )-^2e d ,e t! ( _ -^e d ,o)- 

The 22 topological sectors are given by the con- 
jugacy classes of the pair (91,(72) with the condi- 
tion gi9 2 _1 9r 1 92 = 1, an d are explicitl y 28 ! 29 (1,1), 
(1,-1), (-1,1), (-1,-1), (l,*a a ), (-l,ia a ), (ia a ,l), 
(i(J a , —1), {ia a , icr a ), and (icr a , ~i<J a ), with a = x,y, z. 



D. A Different Perspective 

Previous discussions are based on uniform Af r = z in 
constraint Eq. (jHJ. A different perspective by allowing 
non-uniform Af r will dramatically simplify the picture 
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and results. Do a sublattice-dcpendent orbital rotation 

b A ^{U T -iT y )- x -b A , b B ^(U T -tTy)- 2 -b Bl b c ^b c . 

(52) 

The constraints Eq. (J5J become 

Ma = y, Mb =x, M c = z. (53) 
The mean-field Hamiltonian Eq. (|2T|) becomes 

H\\F = ^2 \f J,n r ~ ' tV(jV ' + ^- c 

3 (54) 
-^2(bl- ir v ® icr y • bl +ed + /i.e.) , 

d=l 

and is a spin-orbital singlet. In this gauge choice the 
quaternion spin liquid state looks like "Tsunetsugu- 
Arikawa orbital ncmatic state" , with T v = 0, T x = 
0, T z = on sublattices A, B, C respectively. This 
"orbital order" reduces the orbital SU(2) gauge free- 
dom to quaternion group. However it does not break 
any physical symmetry in the spin liquid phase. The 
three-sublattice structure becomes physical only upon a 
spin-orbital-entangled condensation of Z at k = 0, with 
(B r ) = Z. The PSG of bosons under this gauge is similar 
to the PSG of Z. In this picture it is clear that the low 
energy theory Eq. (|47p contains the coupling of spinon 
field B ~ (B) to the three Higgs fields Ma,b,c, 

[ d 3 x\x(\TT[B^BT}-Af A \ 2 

J 1 (55) 
+ \Tr[& Bt] -M b \ 2 + \TxpBr\-Mc\ 2 )}- 

Note that B = Z + iW where W is a different SU(2) 
rotor field corresponding to the high energy mode in 
Fig. [2j This gapped mode has to be included because 
Tr[£tgr] = iTr[Z^Wr] - iTr\W^Zr]. 

The structure of the quantum limit wavefunction is 
much more transparent in this gauge. Bond amplitude 
f rr i in the quantum limit becomes 

f , = ( ^ ^ = r±e l,2,3, (5Q) 

Jrr y0, otherwise. 1 ' 

The loop weight Eq. (f5"Tj) becomes 

w t =n l - Tr[T riT yT r2 Ty . . . TrLT y] 

= N L ■ Tr[(-iM* ri ■ r)(-iM* r2 ■ t) . . . (—iAf* L ■ r)]. 

(57) 

(— iAf*-r) = —ir y , —iT x 7 —ir z for r on sublattice A, B, C 
respectively. Then it is very easy to see that the trace can 
only be ±2 or 0. In fact the trace is nonzero only if the 
number of A, B, C sublattice sites on the loop, N a ,b,c[^ 
are of the same parity, N A [i] = N B [P\ = N c [i] mod 2. 



III. DISCUSSIONS 

In Ref. [1^ it was argued that no gauge invariant bi- 
linears of the low energy spinon field can be constructed 
to carry spin-1 quantum number. Therefore it was sug- 
gested the standard projective construction by rewriting 
spin operators into spinon bilinears cannot describe the 
quaternion spin liquid. The argument is indeed true here, 
and there is no gauge invariant spin-1 bilinears of the low 
energy field Z. However there is a high energy branch of 
spinons W (see Fig. [2]) which remains gapped across the 
spin nematic ordering transition [see Eq. (|55[) and related 
discussions]. They together can make gauge invariant 
spin-1 bilinears (Tr[iZ* aW] + c.c - S in Eq. ©). This 
situation was overlooked in the analysis of the low energy 
theory in Ref. l28l . So there is no real contradiction. 

An important issue is what spin-1 Hamiltonian may 
favor this quaternion spin liquid as the ground state. In a 
numerical study of nearest-neighbor bilinear-biquadratic 
Heisenberg Hamiltonian on triangular lattice^!, 

H=J2 JSi ■ Sj + K(Si ■ S 3 ) 2 , (58) 

<ij> 

it was found that a three-sublattice "antiferroquadrupo- 
lar" state, same as the Tsunetsugu-Arikawa proposal^, 
is the ground state if K > J > 0. This may serve as 
the starting point. Farther neighbor and mutiple-spin 
interactions can then be added to destroy the long-range 
order, the sign of these terms may be hinted by looking at 
the loop-products of spinon pairings related to the gauge 
invariant flux^. For example the term — (b i ■ ® ia v ■ 

b 3 ){b) ■ A* k ® ia» • b{)(b k ■ A u ® icr* ■ bt)(b} ■ A* H ®i<?y ■ b\) 
defined on a rhombus ijkl may favor the quaternion spin 
liquid. After projection to the physical spin-1 space it 
contains a term —(Si ■ Sj)(Sk • Si) + (Si ■ Sk)(Si ■ Sj) — 
(Si ■ SiXSj ■ S k ) - (l/2)(Si + Sj +S k + SO 2 similar 
to the 4-site ring exchange of spin-1/2, but with opposite 
sign to that would naturally arise in a Hubbard model^. 

Several possible extensions of the current work exist. 
The two-orbital formalism may also be used to describe 
spin liquids in proximity to other spin nematic or dipole 
orders. One interestin g ex ample would be the Z4 spin 
liquid proposed in Ref. |2!j. A complete classification of 
PSG in this formalism may be a useful guide along this 
direction. 

The projected spin-1/2 Schwinger boson wavefunctions 
have been numerically studied on small lattices by brute- 
force evaluation of pcrmanents^. Generalization to the 
current case is likely very hard, because the overlap be- 
tween a S z basis state and the projected wavefunction is 
not a single but many (2^ ,itB ) permanents. 

The two-orbital AKLT representation can be directly 
generalized to higher spin systems. Spin-S 1 can be repre- 
sented by 2S orbitals of spin-1/2 Schwinger bosons, with 
2S single occupancy constraints generalizing Eq. (|23[) on 
each orbital and S(2S — 1) symmetrization constraints 
generalizing Eq. between each pair of two orbitals. 
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The gauge freedom is [U(1)] 2S x S2S where S2S is the 
symmetric group of degree 2S. This formalism can de- 
scribe higher degree multipole orders by boson conden- 
sation, and spin liquids with even richer gauge structures 
(thus richer topological orders) may be obtained via the 
projective construction. 

Multiple-orbital fcrmionic representation has been con- 
sidered for general SU (N) spins^, and used in the con- 
text of alkaline-earth cold atom systems^. Large- TV gen- 
eralization of the multiple-orbital bosonic representation 
may also be useful in theoretical studies. More recently 
the two-orbital fermionic representation of spin-1 was 
employed^ in hope of describing the experimental evi- 
dence of gapless spin liquid in Ba3NiSb20g^. 

The prototype wavefunction defined by Eqs. (|17H18[ 
I5UII0TI 15M57|) may be of some interest by itself. It is not 
clear how to check directly the quaternion structure with- 



out reference to the mean-field theory. It is also possible 
that confinement happens due to the projection of the 
mean-field wavefunction. The confined phase may have 
nontrivial quantum numbers of the space group^. More 
insight on the amplitude (matrix trace) are much needed 
for these purposes. And it will be very interesting if the 
matrix trace form of the loop amplitudes Eq. (|18[) can be 
used to represent other nontrivial phases. 
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